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Abstract. When considering the statistical properties of a 
bundle of cosmic microwave background (CMB) photons prop- 
agating through space, the effect of 'mixing geodesies' appears 
with a distinct signature that depends on the geometry of 
space. In a Universe with negative curvature this effect is ex- 
pected to produce elongated anisotropy spots on CMB maps. 
We used COBE-DMR data to look for such effect. Based on 
the analysis of a measure of eccentricity of hot spots it ap- 
pears that there is a clear indication of an excess eccentricity 
of hot spots with respect to that expected from noise alone. 
This result must be interpreted with caution as this effect can 
be due in part to galactic emission. If the detected eccentricity 
of anisotropy spots can be attributed to the effect of mixing it 
implies the negative curvature of the Universe and a value of 

a < i. 



1. Introduction 

The detection (Smoot et al. 1992; Bennett et al. 1994) and 
further confirmation (Bennett et al. 1996; Hancock et al. 1994; 
Ganga et al. 1993) of the presence of anisotropies in the cos- 
mic microwave background (CMB), has endowed cosmologists 
with a unique tool for a realistic test of cosmological models. 
The observed large angle AT/T can be produced by a number 
of sources that go from local, to astrophysical, to cosmologi- 
cal effects, among which the most relevant are: emission from 
the Galaxy; gravitational fluctuations produced by voids and 
great attractors; the Doppler effect due to the observer's mo- 
tion relative to the CMB; mixing of geodesic flow; gravitational 
waves; cosmic strings; and gravitational potential fluctuations 
on the last scattering surface - the Sachs- Wolfe effect (for a 
review see White, Scott & Silk 1994). Clearly the complexity 
of determining the cause of CMB anisotropies requires sev- 
eral independent analysis techniques. Aside from the standard 
statistical tests based on the angular auto-correlation function 
and power spectrum, a number of statistics have been proposed 
(Sazhin 1985; Vittorio & Juszkiewicz 1987; Bond & Efstathiou 
1987; Coles 1988; Gott et al. 1990; Martinez-Gonzalez & Cayon 
1992; Martinez-Gonzalez & Sanz 1989) as topological descrip- 
tors to characterize the observed anisotropies: hot spot number 
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density, length and total curvature or genus of iso-temperature 
contours, number of upcrossings, area and eccentricity of hot 
spots, and Euler-Poincare characteristic. The study of these 
topological descriptors is based on the well established theory 
of the geometric properties of the excursion set of random fields 
(Adler 1981). 

The great potential of topological analysis of CMB data 
can be attested by the large number of results obtained thus 
far: the genus descriptor has been successfully used to place 
important restrictions on the shape of the spectrum of pri- 
mordial perturbations P(k) (Smoot et al. 1994; Torres 1994; 
Torres 1995a; Torres et al. 1995); further restrictions on P(k) 
have been obtained from the peak distribution analysis (Fab- 
bri & Torres 1995, 1996); a possible fractal dimension in the 
universe has been detected with the analysis of the contour 
characteristics of anisotropy spots (De Gouveia 1995); studies 
of hot spot number density (Torres 1995b; Cayon & Smoot 

1995) support the results given by the genus analysis; tests of 
Gaussianity of the CMB field have been performed using the 
genus and the correlation function of local maxima (Kogut et 
al. 1996); finally, percolation and cluster analysis (Naselsky & 
Novikov 1995) of CMB maps has been proposed as a powerful 
tool to test the Gaussian character of CMB anisotropies. The 
information carried by the CMB properties can even have a 
direct impact on accelerator physics (Gurzadyan & Margarian 

1996) . 

We propose to use a topological analysis, based on the 
eccentricity or elongation measure of hot spots, to look for 
the signature of the mixing of geodesies effect (Gurzadyan & 
Kocharyan 1991, 1992). Because of the appealing possibility 
of extracting information about the curvature of the universe 
from the latter effect, we will deal with it in some detail. A 
consequence of mixing phenomena in geodesies propagation is 
the appearance of highly distorted anisotropy spots in cosmic 
background maps (Gurzadyan & Kocharyan 1993a, 1993b). In 
the case of an open geometry, anisotropy spots associated with 
a bundle of photons propagating in free space would appear 
as highly elongated hot spots, while in a flat space hot spots 
would show a symmetric circular shape. This visible effect is a 
direct consequence of the geometry of space and as such it is a 
unique tool to discriminate between open and closed cosmolog- 
ical models. However, confusion due to galactic contamination 
at high latitudes and similar elongated patterns produced by 
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instrumental noise and cosmological structure characterized by 
large coherence angles makes this effect very difficult to ob- 
serve. 

Fortunately, at least 3 different phenomena produced by 
the effect of mixing geodesies and associated with the observ- 
able characteristics of CMB, have been predicted: 

(a) Isotropization. The decrease of the amplitude of CMB 
anisotropy after the decoupling epoch: the degree of isotropiza- 
tion, i.e. the fluctuation damping factor can yield Hp to 10 3 for 
n = 0.2-0.4 (Gurzadyan & Kocharyan 1991,1992)0; 

(b) Angular dependence. The behaviour of anisotropy as 
a function of the sky and smoothing (beam) angles, namely, 
the independence of the autocorrelation function on the sky 
angle, and its dependence on the smoothing angle (Gurzadyan 
& Kocharyan 1996); 

(c) Maps. The complex topological structure of the CMB 
sky maps, showing intrinsic threshold independent elongation 
of the shapes of both hot and cold spots (Gurzadyan & 
Kocharyan 1993a, 1993b). 

In this paper we will deal only with the third phenomenon. 
Since, in principle, effects of apparent elongation might occur 
also owing to other physical mechanisms (the same is not ex- 
cluded also for phenomena (a) and (b)), the unambiguous con- 
firmation of the observational discovery of the effect of mixing 
geodesies could be reached by means of a comprehensive study 
of the complete set of predictions. As it is briefly discussed 
at the end of this paper, at least the available evidences do 
not contradict to all these predictions. A confirmation of the 
effect of mixing geodesies would constitute powerful probe of 
the curvature of the Universe. 

2. Effect of mixing of geodesies flow: Maps 

The idea of the effect of instability of trajectories of freely mov- 
ing particles can be most clearly demonstrated via the Jacobi 
equation written for spaces with constant curvature k: 



d n/d X + kn = 0, 



(1) 



describing the behavior of the vector of deviation, n, of 
close geodesies. Solutions of this equation are determined by 
the sign of the curvature: when k < one has exponentially 
deviating geodesies. 

A more rigorous treatment (Gurzadyan & Kocharyan 1991, 
1992) includes the study of the projection of geodesies from 
(3+l)-dimensional Lorentzian space to a 3D Riemannian one, 
and the behavior of time correlation functions for geodesic 
flows on homogeneous isotropic spaces with negative curva- 
ture. Geodesic flows, being Anosov systems (locally if the space 
is not compact), are exponentially unstable systems possess- 
ing: the strongest statistical properties (mixing) , non-zero Lya- 
punov characteristic exponents, and positive Kolmogorov- Sinai 
(KS) entropy h (Arnold 1989). 

For Anosov systems two geodesies in 3-space deviate expo- 
nentially according to the law 



n(A) =n(0) exp(xA), 



(2) 



x In Eq(12) of the interesting paper by Ellis and Tavakol (1994) 
an additional numerical factor 2 is necessary, as readily follows 
from Pesin's theorem; in the rest those authors indeed confirm 
the above mentioned results on the decrease of CMB anisotropy 
via geodesic mixing. 



where x is the Lyapunov exponent. 

For a homogeneous isotropic Friedmannian Universe with 
k = — 1 the Lyapunov exponent is determined by the only 
parameter a, the diameter of the Universe: 



X = l/o, 



(3) 



while Lyapunov exponents vanish when k = 0, +1. 

Time correlation functions, describing the decrease of per- 
turbations also decay exponentially as determined by KS- 
entropy: h = 2\- 

For the Universe expanding as 



a(t) = a(t )(t/t )°', 



(4) 



the relation between the quantitative measurement of the 
distortion of patterns, e, and Q is given by (Gurzadyan & 
Kocharyan 1991, 1992): 



ln(l/e) 

\ i+zn > 



a/(l-a)[l- (1 + Zi) 1 " 1 /"] a< 1 
ln(l + zi) a = l 



(5) 



where Qo is its present value [z = 0), zi corresponds to 
the time when matter becomes non relativistic and z(~ 1000) 
to the decoupling time. We deliberately used the power law 
representation for the expansion law, (instead of trivial use 
of the corresponding exact Einstein equation) to demonstrate 
the role of the expansion rate in the described effect. It is not 
excluded that in the future the geodesic mixing effect can be 
useful to gain information also on the expansion rate of the 
Universe after the decoupling epoch. 

The parameter of elongation e defined via the divergence 
of geodesies in (3+l)-space: 



1(0) a (t ) 

approaches 1 when Q tends to 1 as shown in Fig. 1. 



(6) 



Fig. 1. Elongation parameter e as a function of Q for two exponents 
for the law of the expansion rate of the Universe, a = 1 and a = 2/3. 

The typical pattern of a hot spot as seen today in a k < 1 
universe would exhibit a very complex shape. The elongation, 
e, measures the smallest-to-largest ratio of diameters of one- 
connected regions and is related to the "degree of complexity" 
of anisotropies. 
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3. Analysis of DMR Maps 

COBE's Differential Microwave Radiometers (DMR) have 
mapped the microwave sky at three frequencies: 31, 53, and 90 
GHz (Smoot et al. 1990). The analysis presented here is based 
on DMR's 53 GHz maps (which have the best sensitivity) of the 
four year data set (Bennett et al. 1996), and is a continuation of 
the study based on the one year data-set (Gurzadyan & Torres 
1993). Signal and noise maps were prepared by adding and sub- 
tracting the two independent DMR channels (ie. 0.5A + 0.5B 
and 0.5 A — 0.5B) and Gaussian smoothing (a = 2.9°) the re- 
sulting maps. A Galaxy cut \b\ < 20° has been applied to the 
maps. The geometric characteristics of hot spots are quite sen- 
sitive to galactic cuts below 15° to 20° but beyond 20° our 
results are stable. 

DMR maps are digitized on a grid of 6144 approximately 
equal area pixels of size ~ 2.9° (Torres et al. 1989). A hot 
spot is defined as a continuous region of the map formed by 
pixels whose temperature is higher than a preset threshold 
T v = vtj, where a is the standard deviation of the sky tem- 
perature (monopole and dipole subtracted). The algorithm to 
identify hot spots (Torres 1994, Torres 1995a) relies on the con- 
struction of binary tree structures out of the set of connected 
'hot' pixels (for a given temperature threshold). The number 
of pixels in a tree gives the area of the spot. The eccentricity 
parameter e l v = r\ ain lr % max of the i-th hot spot at threshold v 
is found as follows: first, the center of the spot is estimated as 
the point with x c , y c coordinates equal to the mean of all the x 
and y coordinates of the pixels that form the spot respectively; 
second, the distance from this center point to all the pixels 
that lay on the boundary of the spot are calculated; r % max and 
r min are t ne largest and shortest of these distances. The above 
mentioned procedure is repeated for several threshold levels 
(y : 1.0 — 3.0 in steps of Ai/ = 0.25). The geometric descrip- 
tor e„ used to measure the elongation of anisotropy spots at 
threshold v is the average of all the ej, at a fixed v. 

At this point one can question whether the defined e v re- 
ally corresponds to the theoretical parameter of elongation 
in Fig. 1. The correspondence of the Lyapunov numbers or 
KS-entropy with parameters observed either with experiments 
or computer simulations is an essential problem. However, it 
has been empirically established that a correspondence with 
macroscopic parameters does exist (at least qualitative) and 
just this fact determines the role of those parameters as im- 
portant tools for the study of nonlinear phenomena (Hilborn 
1994 and references therein). Another question is whether the 
specific procedure of estimation of the mean elongation of spots 
can be crucial or not. Obviously other algorithms in principle 
can lead to numerically not absolutely equivalent results, but 
a drastic change of the situation seems unlikely, especially at 
high thresholds (Sommerer & Ott 1993 and references therein). 

Figure 2 shows e v for the sum and difference DMR maps 
for threshold levels in the range v : 1.0 — 3.0. Data for v < 
1.0 is dominated by noise and for v > 3.0 is limited by the 
large statistical error due to the small number of spots at high 
thresholds. 

4. Results 

The signature of the mixing of geodesies effect is a clear 
one: hot spots should have a fixed eccentricity independent 
of threshold level, and if these hot spots are elongated it is an 
indication of mixing in open spaces. In principle it is possible 



Fig. 2. Measured eccentricity parameter e„ for the COBE sum maps 
(solid) and difference maps (dash) compared with the expected ec- 
centricity for noise Monte Carlo realizations (points with error bars) . 

to detect eccentricities smaller than 1, the problem however 
is to disentangle the effect due to mixing from instrumental 
noise which by its stochastic nature is expected to produce 
anisotropy spots with certain degree of elongation. 

In order to evaluate the statistical significance of a possible 
detection of the mixing of geodesies effect we have performed 
Monte Carlo studies of noise maps that take into account in- 
strumental noise, COBE's sky coverage and pixelization. Noise 
maps were generated by assigning to each pixel on the map a 
temperature equal to a random number extracted from a Gaus- 
sian distribution with dispersion a = a\j\fN where <ti is the 
corresponding sensitivity for one observation and N the num- 
ber of observations. 

Maps for both A and B channels were generated indepen- 
dently. A difference map was formed and Gaussian smoothed 
just as it is done with the real data. The same algorithm used 
to obtain the eccentricity parameter of the COBE maps was 
used for each one of the Monte Carlo noise realizations. Fig- 
ure 2 shows the Monte Carlo mean eccentricity and the 1-a 
error bars expected from noise maps. Knowing the eccentric- 
ity expected from noise, one can estimate the probability that 
the observed eccentricity parameter can be produced by noise 
alone. Table 1 gives e„ for the COBE sum and difference maps, 
and the deviations from the mean eccentricity of Monte Carlo 
noise simulations. 

The x 2 statistic computed with the 9 data points in the 
range v — 1.0 — 3.0 and the corresponding noise Monte Carlo 
points is 5.6 and 73.0 for the (A — B) and (A + B) maps respec- 
tively. The low x 2 associated with the difference maps was ex- 
pected and is an indication of the accuracy of the Monte Carlo 
simulations. On the other hand, the high \ 2 obtained when 
the data from the signal maps is compared with noise data is a 
clear indication of an actual detection of elongated anisotropy 
spots. The measured eccentricity parameter at threshold levels 
1.5 and 2.5 in particular exhibit a > 4-a deviation with respect 
to noise. The average deviation in terms of standard deviations 
of e„ +B from the corresponding Monte Carlo result for the 9 
bins considered here is 3a. From the mixing of geodesies ef- 
fect one would expect a constant eccentricity for all threshold 
levels. Due to the large dispersion of the measured e v values 
it is not possible to make a strong statement in favor of the 
hypothesis of a constant e„ independent of v. However, it is 
seen that the contribution to the \ 2 is roughly the same inde- 
pendent of v. Under the hypothesis of a positive detection of 
elongated anisotropy spots and using the 9 points in Table 1, 
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Table 1. Eccentricity parameter of hot spots on COBE maps 
(e A+B , e^~ s ) and comparison with Monte Carlo noise maps 
(e™ c ). A A+B and A A ~ S denote the difference (in standard 
deviations) between the measured eccentricities and the mean 
eccentricity of noise Monte Carlo maps. 



V 


,A+S 


A-B 


MC 


A A+B 


A A-B 


1.00 


0.480 


0.570 


0.582 


3.551 


0.447 


1.25 


0.547 


0.590 


0.617 


2.271 


0.886 


1.50 


0.497 


0.616 


0.646 


4.288 


0.881 


1.75 


0.576 


0.691 


0.674 


2.432 


-0.429 


2.00 


0.600 


0.686 


0.699 


2.049 


0.267 


2.25 


0.595 


0.819 


0.719 


2.060 


-1.675 


2.50 


0.427 


0.685 


0.741 


4.017 


0.719 


2.75 


0.574 


0.732 


0.761 


1.752 


0.272 


3.00 


0.507 


0.705 


0.776 


1.815 


0.480 



the measured eccentricity is e = 0.53 ± 0.06. 
5. Discussion 

An excess elongation as possible genuine feature of the hot 
spots on COBE maps is well established at least at threshold 
levels v between 1.0 and 3.0. The detected excess elongation 
shows a statistically strong independence on the threshold, con- 
trary to the case for the difference (A— B) maps where e„ shows 
a clear correlation with threshold. This characteristic depen- 
dence of the eccentricity parameter with threshold for noise 
maps is clearly manifest in Monte Carlo noise realizations. If 
the detection of high eccentricity of hot spots is attributed 
to the effect of mixing of geodesies flow, our analysis implies 
fl < 1. However one must also consider the fact that the de- 
tected effect can be produced in part by galactic contamination 
of the maps. 

The numerical value of the obtained elongation parame- 
ter might have been obviously affected by a number of effects 
of cosmological and non-cosmological nature. Primordial fluc- 
tuations on the surface of last scattering, for example, would 
also produce non-zero ellipticity even in flat universe (Bond 
& Efstathiou 1987). However, due to the stochastic nature of 
this effect, the dependence of the ellipticity parameter with 
threshold would be just like that of noise, contrary to what we 
have observed (i.e. a threshold independent eccentricity). More 
statistics and higher accuracy is necessary for a reliable numer- 
ical evaluation of 51. A question thus arises as to how one can 
distinguish whether the observed elongation is really due to 
the effect of geodesies mixing? Indeed, although the "geodesies 
mixing" has put forward the idea of looking for the threshold 
independent elongation of hot spots on CMB sky maps as a 
signature for Q < 1, there are other mentioned effects that 
can be considered as well. Thus, to answer the question posed 
above it is important to have a thorough consideration of other 
observable consequences of the effect of mixing geodesies. As it 
is mentioned above, among those consequences is the decrease 
of the amplitude of initial fluctuations during the expansion of 



the Universe; the expected decrease depends on the curvature, 
the value of Q and the time of last scattering (Gurzadyan & 
Kocharyan 1992). 

More interesting, however, can be the property concerning 
the dependence of the temperature autocorrelation function 
C(9, j3) on the sky angle 6 with a beam size (3. It was shown 
(Gurzadyan & Kocharyan 1996), that in negatively curved 
spaces C(8,/3) tends to become constant with respect to 9 (i.e. 
independent on the sky angle) independent of the initial per- 
turbation spectrum, but should depend on the beam size (5 by 
the asymptotic law 

C(9,f3) _ const 
f 2 P ' 

In other words, the observation with different beam sizes 
should lead to different values for the amplitude of the CMB 
anisotropy, namely the measured anisotropy decreases with in- 
creasing beam size. 

Regarding the latter effect, note that some measurements 
by smaller beam sizes (Netterfield et al 1995; Ruhl et al 1995) 
seem to indicate higher values for the anisotropy amplitude as 
compared with the COBE data, even though one has to dis- 
tinguish the contributions of other effects such as the Doppler 
peak for example. 

The peculiar feature of the effect of geodesic is its crucial 
dependence on the curvature of the Universe, i.e. its disappear- 
ance in flat and positively curved spaces, and independence on 
various models of dark matter or the spectrum of initial per- 
turbations, etc. 

Obviously nothing is happening to an individual photon 
during the free propagation after the last scattering epoch, 
and these effects are entirely statistical and determined by the 
principal limitations of obtaining information during the mea- 
surements, i.e. by the impossibility of the reconstruction of the 
trajectory of the photon while observing within finite smooth- 
ing angle and time period, due to the overlapping of exponen- 
tially deviating geodesies in any given cut of phase space (for 
detailed discussion of this and similar problems see Gurzadyan 
& Kocharyan 1994). Note that these effects differ from those 
of classical chaotic Mixmaster models, and are alternatives to 
inflationary scenarios (Gurzadyan & Kocharyan 1994; Cornish 
et al 1996). The measured effect of elongation of anisotropy 
spots in CMB sky maps can be the direct footprint of the neg- 
ative curvature of the Universe and therefore has a considerable 
impact on cosmology. 
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